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HEAT PIPE STABILITY I: A PRELIMINARY INVESTIGATION 
INTO THERMALLY ASSISTED CAVITATION 


ABSTRACT 

This paper introduces in a preliminary manner the notion of 
thermally assisted cavitation by localized fluctuations in capillary 
forces. Because cavitation in liquids can be closely approximated 
by an isothermal process, only momentum and mass balances are 
used to introduce the notion for liquids under mechanical tension. 
The corresponding phenomena for the electrocapillary effect is not 
considered. Brief attention is given to developing a stability theory 
in terms of the stiffness and compliance coefficients for a working 
fluid. Interestingly, the particular thermodynamic approach taken 
can be used to suggest experiments relating working fluid perform- 
ance to meniscus behavior. 

A correction to the liquid flow equation is suggested. 
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HEAT PIPE STABILITY I: A PRELIMINARY INVESTIGATION 
INTO THERMALLY ASSISTED CAVITATION 


INTRODUCTION 

In referring to the concept of wetting, it is customary to relate the referral for 
a given material to a contact angle measured at room temperature. | The angle 
so measured is therefore an equilibrium value; variations in thermodynamic 
forces will cause variations in this angle, but the fluctuation normally disappears 
when standard conditions are reestablished. Engineering calculations custom- 
arily assume a contact angle of zero in the Laplace equation because angles are 
usually judged small. If this assumption is justified in all cases, then the equi- 
librium state of the meniscus is defined for zero angle, and any departure result- 
ing from fluctuations in thermodynamic forces causes a shift away from equilib- 
rium with the possibility of an accompanying shift in contact angle away from 
zero. A contact angle of zero means simply that the adhesive forces between 
liquid and substrate exceeds the cohesive forces between liquid particles suffi- 
ciently to cause spreading, and the liquid to flow sufficiently to form a layer 
approaching the thickness of a monolayer, that is, the sine of the contact angle 
is zero, and the covering layer at prevailing conditions appears to the substrate 
as an extension itself. 

Because the Laplace equation forms the basis for steady state momentum and 
energy transfer generated by capillary phenomena, it must also serve as a starting 
point for stability investigations. We commence our discussion from the first 
variation, viz, 

S'V 5 J? y 

S(P - P, ) = — j- cos 9 7 cos 6 — ^ sin 6 6 6 > ( 1 ) 

® K K 

since ultimately it is the functional relationship of pressures, radii of curvature 
etc. with space and time coordinates, which establishes the notion <|>f stability 
or lack thereof. Ideally, it is possible to transfer thermal energy isothermally 
because density fluctuations across an interface at the phase change temperature 
generates a pressure fluctuation and fluid flow. However, because of viscous 
effects, fluctuation dissipation occurs and interferes with energy transport. As 
a result, temperature gradients are normally introduced, and operation is re- 
moved from equilibrium considerations becoming subject to analytical for- 
mulation as a steady state problem. Locally, however, to speak of processes 
as occurring isothermally offers convenience, and is justified by treating a de- 
vice such as a heat pipe as being made up of sections, and with some lack of 
consistency by saying also that thermal gradients are small. A gradient may be 
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small, but if the dependent variable is not, linearization may not be justified, 
particularly if the structure of stability relationships is required. After this 
ado, we begin by examining equation (1) for the isothermal case. Since 88 <0 
for 0 0, then 5(P g - P L ) = O provided the curvature compensates such that 

5R > O, or at the boundary the meniscus becomes tangent. Changes have there- 
fore occurred in a conservative manner. Keeping this conservation relationship, 
the isothermal requirement may be relaxed, 8 ( P g - P L ) = 8y /R (cos0), and over 
some temperature range, where 8y is significant, pressure fluctuations are con- 
trolled by surface tension whereas 5R is controlled by adhesion to the wall. Next, 
looking at the situation of fluctuations about 0 = O, 5(P g -P L ) = O and the func- 
tional relationships are conservative if along some segment not attached to the 
boundary, the meniscus changes shape such that they are related by 5R > O then 
Sy > O and AT < O, or the meniscus complies by becoming more rigid by re- 
moving particles from the liquid. Reversing equality signs produces the oppo- 
site effect. In the isothermal case, pressure changes must occur for 5R. For 
5R > O, 5 P l > 5P g and the density of gas over the given locale in the meniscus 
will increase; for SR < O condensation occurs. As a basis for the above discus- 
sion, pinning of the meniscus at the boundary at 8 = O was assumed, and there- 
fore a film of liquid, a monolayer or so in thickness, exists along some portion 
of the boundary. Also at this contact angle, a gradient in heat of vaporization 
exists from the film to the meniscus proper. 

Although the above discussion of the Laplace equation is incomplete, enough has 
been said to provide some basis for the elaboration of the stability problem and 
its relationship to the bubble formation mechanisms outlined below. 


Approach 


The steady state transfer of heat by a liquid under tension occurs because a 
differential in curvature exists between evaporation and condensation interfaces. 
To zeroth order, the difference is a design feature incorporated into a device 
such as a heat pipe. Flow velocity is as a result dependent upon a pressure 
gradient defined relative to the distance between heat source and sink. A per- 
formance description of say a heat pipe can be adequately represented to this 
approximation provided the stiffness and compliance coefficients governing be- 
havior are able to react to displacements from steady state in a manner such 
that the thermodynamic forces governing behavior are maintained at their mean 
value; that is, a fluctuation must decay as exp(-At). If this is not the case, per- 
formance is determined by local conditions , particularly the states of local 
meniscii. Now, since viscous effects — both transverse and longitudinal — are 
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responsible for controlling the dissipation of mechanical perturbations at any 
temperature, it is appropriate that the momentum balance equation serve as a 
starting point, viz, 


with 3 t and, j serving to specify that a time derivative of the i-th velocity com- 
ponent is to be taken, and a coordinate derivative of the pressure tensor Py ; 

Fj represents a component of body force. Inertial terms of form pvj v^ j have 
been neglected because motion is presumed slow, and between boundaries closely 
spaced. Taking the variational derivative, multiplying by the velocity variation 
and integrating over volume elements dV, the result is 


i»i«V 2 = - «Vi sPjj 


J*\ i sp ij 


dV. 


Incorporating the definition of the pressure tensor 


p ij = p 5 ij + " ,j. 


that is, its decomposition into static and viscous parts, the result is 


( 2 ) 


(3) 

(4) 



+ 


(5vj) 2 dV = | fiVjSp, + fiv i -2.6(v ij j + vj ii ) | 
f 5v i,i 5p i dV-ir ? J'(5v jJ + 6v j>j ) 2 dV. 


(5) 


Under steady state conditions 5 v ; = O at the end points, and <5 vj s = O because 
the fluid is incompressible at the temperatures in question, and the fluctuation in 
kinetic energy of the fluid is always damped out by viscous forces: 


(5vj) 2 dV = J*(5v t j + 5vj j) 2 dV. 


In the case of tension controlled fluid transport at temperatures where viscosity, 
surface tension, heat of vaporization, and vapor pressure are varying functions 
of temperature, Svj at the meniscii in question need not be zero, and 5vj ; may 
as a result differ from zero. Under these conditions, stability related to the 
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behavior of all liquid-vapor-inert gas interfaces in question may effect the bulk 
liquid in the vicinity of a portion of these interfaces. Considering terms in (5), 
other than the Rayleigh term, stability or the damping of velocity fluctuations 
is not apparent unless each of the four terms is considered separately. The 
first two terms may be written in explicit form: 


5Vj • 5p t 


=2 /svj6pj 

jKi J 


dAj K and 


Sv l2 <(T i .) +V i ■i ) | i / 5V i 2 S(V i.J +, j.i ) 


dA: 


iK' 


( 6 ) 

(7) 


Rewriting equation (5) after rearranging such that the terms responsible for 
instability are to the left of the equals sign and the Rayleigh term to the right 
of it reads 

fi d t( 5v i) 2 dV -/ 5v i • 5 Pi • d V -/sv^StVi j+Vj jldA^ 

r j f ^ 

-JSvj i • 5pj dV - -~n J (5v i;j + 5v jti ) 2 dV. 

The first term is controlled by viscous forces whenever the change of acceleration 
with time is negative; however, instability sets in if both 5vj and 9 t 5vj are nega- 
tive, that is the flowing mode of the system ceases to function, hi referring to 
the following two terms, it is to be remembered that not only is total area a 
factor, lets say the areas Aj k across which hydrostatic pressures are exerted, 
but also areas A iK which contain the flow, but also local areal fluctuations. In 
general for any functional device such as a heat pipe, interfacial area within 
theisystem is not a conserved quantity, although total volume is. Observed from 
this viewpoint, a tension controlled device is a metastable one at best, excepting 
perhaps a system using a liquid metal where the temperature difference between 
melting and boiling points is large. Clearly if the aggregated sum of the varia- 
tional pressure terms in (8) exceeds the Rayleigh term the system is unstable. 

Looking into the first of these terms, dA^, 5Pj and 5Vj must all be > O, if per- 
formance is to be increased. For meniscii neighboring a heat source, the 
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aggregated increase in areas must be small and related to 5pj» which 
may be defined, e.g. , 


5 Pi 


a y cry a 

V*; 1 


5R: 

R i 2 


sin 6, 50:. 

K i 


(9) 


the factor, 5v if simply represents a functional velocity variation for the Spj 
given in (9). For the second term, the shear term may be written in the form 




dy 

dT 


AT 


( 10 ) 


dy 

1 dT. 1 

dT 

O 

Cr 

J* 

1 


dy 

dT 


dT dxi 

„ E, as is ’ 


where variations in temperature gradient have been considered relative to the 
direction of flow, SXj , and along an element of arc, 5 S. Both terms in general 
are greater than zero: 


e| 5 v s Pi 

jKi ' 


dAj K > O, 


(ID 


£ f Sv i 1 ( V i, j + V j, i> (1A iK >0 


( 12 ) 


and are destabilizing. This observation may be seen in another manner by con- 
sidering the Gibbs free energy, G. Because 5pj > O, there must be an increase 
in flow work and therefore 



5G > 0, 

(13) 

that is 

6G = v5p + 7 dA, 

(14) 

or 

N5p = v5p 



because only isothermal variations are being considered although flow' may be 
taking place above some reference temperature level. But since flow 7 is 
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continuing 'with 50 < O, then 


d^>0 (15) 

and the available energy per particle transferred is increased. 

Turning now to a discussion of the third term or dilational factor, we commence 
from the equation of continuity. At the surface of evaporation, the continuity 
condition may be written 


= P L Au, (16) 

n v 

with QA V representing rate of mass transfer based on thermal properties, and 
p L Au mass transfer rate relative to the fluid properties. Based on the Laplace 
equation, the surface shape is fixed; once, however, a generalized force such 
as a thermodynamic force is applied in the neighborhood of a surface, the sur- 
face becomes only quasi-fixed, and as a result fits a free surface condition 
with the initial shape given by the Laplace equation, 

V = v 0 + v } , (IV) 

but rj = rj (y, t) or 


from (17) 


dr? dr} 

u = Tt + v ^ 


9t?, 

u = ^T +V 


3t? 0 3r? t 
9y + dy 


(IB) 


(19) 


that is, the particles composing the surface must have a velocity component 
normal to the flow direction which may or may not coincide with meniscus 
velocity in that direction. From (16) 


^Lv(Mu) = C 

but since 5v 4 are not zero for all interfaces, and points in the bulk liquid, we 
may write for a surface from (20) 


if 


P 2 A 2 u 2 p^AjUj 


3(pA) 

at 



o, 


( 21 ) 
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and a discontinuity in flow has appeared in the neighborhood of the surface. For 
simplicity write 


or 



If the bulk liquid does not comply then 


( 22 ) 


A (^-) = P g 4rr R2 R (23) 

9 

with u - p 2 + R and a bubble forms. Under circumstances where bulk liquid 
complies then the equation of continuity for bulk, 


or 

and 


~^ + V(pAu) = O 

(24) 

pA(u 2 -u,) ^ O, 

(25) 


(26) 


implying that there is an adjustment in area such as a bubble at some point in the 
bulk liquid not in the immediate vicinity of the meniscus. Whether or not a 
tension- simulated bubble so formed is a cavitation phenomenon resides in whether 
or not the process is isothermal. This aspect may be checked from an energy 
balance in the neighborhood of the bubble: 

p L C Pl 4ttR 2 A $ AT - h v p g |?rR 3 (27) 

and 

A£ = [DAt] % (28) 

with D = (K/ pCp L ) and p g the gas density in the bubble. For water at 15° C the 
temperature drop amounts to less than a degree (1). The above equation states 
that the temperature rise due to thermal diffusion from the interface to a dis- 
tance A£ in At is in this approximation counterbalanced by the vaporization 
of liquid to a radius R, If AT is small then the bubble is produced by tension 
through competing processes resulting in a. small AT and cavitation. On the other 
hand, if AT is not small either in a positive or negative direction, the bubble has 


7 



been formed by nucleation. In either case the left hand side of (8) may be rewritten, 


J d t ( 5 Vi ) 2 d V = \p f 3 t (5R) 2 dV, 


(29) 


with the volume of integration being taken from the surface of the bubble to all 
neighboring meniscii for which 5p i ¥= O. 

Ultimately, of course, stability must be related to stiffness coefficients for the 
working fluid; that is, the derivatives of extensive quantities with respect to in- 
tensive ones. Derivatives of intensive quantities with respect to extensive ones 
serve to define compliance coefficients. The two sets of quantities do not bear 
a simple reciprocal relationship because different variables are held constant 
in each definition. The approach taken below is based upon the requirement that 
the Gibbs free energy be at a minimum. 


For the gas and liquid phases to be at minimum free energy, the second varia- 
tion must satisfy the inequality 


with 


5 2 G > O 


(30) 


5 2 G - 2 
». j 



A£j 


(31) 


and = T, P, N, A, the temperature, pressure, number of particles in the 
liquid phase, and some interfacial area subject to variation. The conservation 
relations 


V g + V L = V (32) 

and 

N g + N L - N (33) 

hold where g and L refer to gas and liquid phases, respectively. On the other 
hand interfacial areas are not conserved, 

° < S V < A m a x- (34) 

Next (31) is diagonalized by completing the square (2) to give 

5 2 G = S Xj if > O. (35) 
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the fi do not have physical meaning, but for (35) to be positive definite each A, 
must be greater than zero for all arrangements of terms used as a starting 
point in completing the square. Results given below are obtained for the order- 
ing: (T, P, N, A). The first two coefficients are 



X 3 and A 4 are given in the appendix because of their complexity, and not described 
in detail here. Equation (36) is clearly > O and uninteresting. Relationship 
(37) on the other hand is interesting so long as 



(38) 


The system is provisionally stable at best, since all possible A's must be greater 
than zero for absolute stability. Letting 


P 


ay 
"R 5 


(39) 


(-) 

Wtna 


(£)(!?) I 


TN A 


so that 



(40) 


(41) 


where V = V L , liquid volume because / V g dP g ~ O. To gain some insight 
into the meaning of this factor, reference is made to Figure 1, which shows the 
volumes and curvatures in question. 
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Figure 1. Illustration For X, 


Referring to Figure 1 , 


av l > o, 

dR < O 

giving 

( 3 ,„. < • 

as required for stability. Stability in this case meaning continued flow even 
though a bubble has formed to conserve area and give momentarily an increase 
in flow. We are now in a position to reverse the description and explain what 
has really happened: A temperature has been reached where the meniscus (<p 0 , 
S 0 ) has cavitated liquid nearby because 



or the vapor pressure (and heat of vaporization) of the liquid is changing more 
rapidly with temperature than surface tension; hence the meniscus cavitates the 



liquid with d - constant in an effort to maintain flow continuity at the liquid 
vapor interface. What appears as increased performance by bubble-flow sta- 
bilization, may eventually, therefore, lead to a break down in flow when the 
bubble grows to sufficient size to cause the meniscus to revert from state (1) 
to (0), making 



(43) 


and the system unstable. 


To lend insight into the order of magnitude of (dV/9P) TN A , 
is shown for liquid ammonia, writing 



■ ; (!) 


N A 


a simple calculation 


and 


/9S\ 1 3h v 1 , AT , 1 , .T 

(ii) PNA '717-7A AT ~^ AT ' 


then at 60° F 


if^v = i.i 

p 4 \ 9T/ 


= 1,17 X 10 -9 (ft 6 /(lb. mass) z (°F) Z )> and 


2 f°-c\ 2i 


relative to 0° F 


-1 h AT = 1.1 X l(T l BTU/(lb. mass)(°R) 2 ) 
T 2 V 


giving 


/ d( 1 /p)\ 

\ ^ / TNA 


> 1.06 X 10' 8 


for stability. 


Finally, for situations where a thermal gradient exists between a heat sink 
and heat source connected by say a fluid artery, the existence of this gradient 
may serve to alter the boundary conditions determining the flow pattern. Nor- 
mally, say, for one dimension, 



d 2 u 
dy 2 


(44) 
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with 



£-0, y = 0 

(45) 

and 

u = 0 at the boundaries. 

(46) 


On the other hand, for a tension controlled boundary, the integration constants 
are obtained from 



(47) 


at the boundary, and / pudA which determines the mass flow rate. Under these 
conditions it turns out that 



at the axis, but it may not be used as a condition in view of (47). 


Remarks 

In a heat pipe, lets say, many interfaces exist: liquid- vapor, liquid-metal, 
vapor-metal. The liquid- vapor area, in particular, may be characterized by 
a variety of curvatures. The fact that the liquid-vapor interfacial area is not 
a conserved, although extended, quantity means that for some range of tempera- 
ture the thermodynamic and mechanical properties of the liquid can cooperatively 
remove it from a state of minimum Gibbs free energy, only to be returned by a 
fluctuation in liquid-vapor interfacial area. 

The approach taken in relating the minimization of the free energy to liquid prop- 
erties does not give a unique solution because diagonalization by completing 
the square of the polynomial can be carried out in twenty-four ways for the phys- 
ical situation considered here. Only one such arrangement of terms is considered. 
However, as can be seen from the discussion of \ 2 , at least a portion of the A 
factors can be used to suggest experiments for investigating meniscus behavior. 
The total number of A's comprises a set of ninety-six. Their determination 
however, may be easily carried out. 

It has been shown also how local meniscus behavior can effect vacinal bubble 
formation, in particular, formation by cavitation, and how such an event is 
effected by viscosity. 
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Extending the notion of stability to rate of entropy production was not undertaken. 
This approach, however, could be more rigorous. Similarly, looking into the 
rate of approach of the Gibbs free energy to equilibrium would be useful, since 
knowing whether or not a fluctuation is damped, or amplified, and under what 
conditions is important, although not critical for an insight into the motion of 
stability. 
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